Dephasing due to background charge fluctuations 
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In quantum computation, quantum coherence must be maintained during gate operation. How- 
ever, in physical implementations, various couplings with the environment are unavoidable and can 
lead to a dephasing of a quantum bit (qubit). The background charge fluctuations are an important 
dephasing process, especially in a charge qubit system. We examined the dephasing rate of a qubit 
due to random telegraph noise. Solving stochastic differential equations, we obtained the dephasing 
rate of a qubit constructed of a coupled-dot system; we applied our results to the charge Joseph- 
son qubit system. We examined the dephasing rates due to two types of couplings between the 
coupled-dot system and the background charge, namely, fluctuation in the tunnel coupling constant 
and fluctuation in the asymmetric bias. For a strong coupling condition, the dephasing rate was 
inversely proportional to the time constant of the telegraph noise. When there is fluctuation in 
the tunnel coupling constant, Gaussian decay occurs in the initial regime. We also examined the 
rate of dephasing due to many impurity sites. For a weak coupling condition with fluctuation in 
the asymmetric bias, the obtained dephasing rate coincided with that obtained by the perturbation 
method using the spectral weight of a boson thermal bath, which is proportional to the inverse of 
the frequency. 

PACS numbers: 03.65.Yz, 73. 21. La, 73.23-b 



I. INTRODUCTION 



Efforts to implement quantum computation have re- 
cently intensified. The application of the quantum bit 
(qubit) to solid-state materials, such as superconducting 
Josephson junctions [1] and quantum dots [2-4], is par- 
ticularly promising, because these implementations have 
the advantage of scalability. In a coupled-dot system, for 
example, the localized states in the left and right dots 
are treated as a basic two-level system, in which the tun- 
nel coupling between the two dots constructs a quantum 
superposition of the dots. This superposition manifests 
itself in coherent quantum oscillation (Rabi oscillation), 
and transition can be induced between the superposed 
states [5]. A nanometer-scale superconducting electrode 
connected to a reservoir via a Josephson junction consti- 
tutes an artificial two-level system in which the charge 
states, coupled by tunneling, differ by 2e, where e is the 
electronic charge. This system has shown clear Rabi os- 
cillation [1]. 

Quantum coherence must be maintained during quan- 
tum gate operation. Dephasing, characterized by the de- 
phasing time, originates from various couplings between 
the qubit and the environment. When the qubit is im- 
plemented in a solid-state system, the effects of phonons 
and electromagnetic and background charge fluctuations 
(BCFs) are important in the dephsing process. The effect 
of phonons has been examined in semiconductor quan- 
tum dots as the source of the dephasing accompanying 
dissipation [6] ; the effect of electromagnetic fluctuation in 



Josephson junction qubits has been extensively studied 
[7]. However, BCFs have not yet been examined sys- 
tematically, in spite of their importance in the dephasing 
process. 

BCFs have been observed in many systems [8-11]. In 
nanoscale systems, they are the electrostatic potential 
fluctuations due to the dynamics of electrons or holes 
trapped at impurity sites. In particular, the charge of a 
single impurity fluctuates with the Lorentzian spectrum 
form, which is called "random telegraph noise" in the 
time domain [10,12]. The random distribution of the po- 
sitions of such impurities and their time constants lead to 
BCFs or 1/f noise [13]. In solid-state charge qubits, these 
BCFs lead to a dynamical electrostatic disturbance and 
hence the dephasing. The effect of 1 /f noise on a charge 
Josephson qubit has been examined theoretically - the 
interaction between the qubit and environment has been 
treated by the perturbation method [14,15], by Gaus- 
sian approximation [16] and by the path integral method 
within a spin-Boson model [14,15]. When fluctuating im- 
purities exist in the substrate, not in the junction [11], 
the coupling between the qubit and BCFs is weak, and 
the perturbation method is sufficient. However, when the 
interaction between the qubit and environment is strong, 
methods that go beyond perturbation are needed. 

In this study, we investigated how the electrostatic dis- 
turbance of time constant To coming from a single impu- 
rity affects the quantum coherence of a qubit irrespective 
of the strength of the qubit-impurity coupling. We also 
examine the effect of many impurity sites. This approach 
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is in clear contrast with previous ones [16,15], in which 
the phenomenological spectral weight of the boson ther- 
mal bath was used to characterize the effect of BCFs. 

We consider two types of couplings between the qubit 
and the environment: pure dephasing and dephasing 
accompanied with relaxation of the population [17]. 
In symmetrical couplcd-dot systems, the former corre- 
sponds to a fluctuation in the tunnel coupling constant, 
and the latter to that in the asymmetric bias [18,19]. 
The mapping from a couplcd-dot system to the Joseph- 
son charge qubit is discussed in Sec. VI. By using the 
method of stochastic differential equations, we obtain an- 
alytically the dephasing rate, which is shown to be always 
smaller than Tq 1 . It should be noted that this dephas- 
ing process does not mean the qubit becomes entangled 
with the environment, but rather it means the stochas- 
tical evolution of an external classical field, suppressing 
the off-diagonal density matrix elements of the qubit af- 
ter being averaged over statistically distributed samples. 

Section II defines the Hamiltonian of the system. Sec- 
tion III explains the method of stochastic differential 
equations. The fluctuations in tunnel coupling and asym- 
metric bias are examined in Sees. IV and V, respectively. 
Section VI is devoted to discussion, including the effect 
of many impurities. Section VII summarizes the paper. 



II. HAMILTONIAN 

The qubit and the effect of a single impurity are ex- 
amined in terms of the following Hamiltonian: 



H = 

H q b = 



Hqb ~\~ Hqb—inipj 

hA 



he, 



(c[c R + c R c L ) + —{c\c L - c R c R ), 



(1) 
(2) 



where 4 R and cl,r are the creation and annihilation 
operators of the left and the right dots, assuming a sin- 
gle level for each dot, as shown in Fig. 1(a). The A is the 
tunnel coupling between the dots, and e is the asymmet- 
ric bias between them. The interaction between a qubit 
and the charge at the impurity site is described by the 
following Hamiltonian: 



H, 



qb—imp 



Ti.J T 



+ ^-(c{c L -c R c R )2(dU-l/2), (3) 

where Jt is the magnitude of the fluctuation in the tun- 
nel coupling, Jb is the magnitude of the fluctuation in 
the asymmetrical bias, and d) and d are the creation and 
annihilation operators of the charge at the impurity site, 
respectively. 



(c{c R + c R c L )2(dU-l/2) 



qubit 



B 



(a) 
X(t) 




1 

o 
-1 



(b) 



FIG. 1. (a) Schematic diagram of coupled dots and en- 
vironment constituting an impurity site and the electron 
reservoir, (b) Example time sequence of random telegraph 



We assume that the time evolution of statistical vari- 
able X(t)(= 2 < Sd — 1/2 >) is a Poisson process. As- 
suming a strong coupling between the charge at the impu- 
rity site and the nearby electron reservoir, the dynamics 
of the charge induces not Gaussian white noise, but ran- 
dom telegraph noise [12,18,19]. We therefore consider the 
effect of random telegraph noise with characteristic time 
constant r , where the statistical variable X{t) takes the 
value 1 or -1 (Fig. 1(b)). 

The time constant is determined by the barrier height 
of the electron or hole trap and the temperature, like t = 

Ae W/k B T^ where w ^ fcs; T> and A 

arc the activation en- 
ergy of the impurity potential, Boltzmann's constant, the 
temperature, and the temperature-independent prefac- 
tor, respectively [20]. When the temperature decreases, 
To becomes longer. The telegraph noise has been exper- 
imentally observed to take a value of 1 or -1 with asym- 
metric probabilities, which arises from the difference be- 
tween the Fermi energy of the electron reservoir and the 
energy level of the impurity sites. [12]. To include this 
asymmetric weight, we introduce asymmetric probabil- 
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ities p u and pd, which relate the asymmetric transition 
rates for the process from -1 to 1 ( t^ 1 = (p u To) _1 ) to 
those of the opposite process ( r^ 1 = (p^o) -1 )- 

We neglect the backaction from the qubit to the charge 
at the impurity site, so this environment reduces to a clas- 
sical stochastic external field. We also assume that the 
temperature is high enough for the effect of the quantum 
fluctuation of the charge between the impurity site and 
electron reservoir to be neglected. 

For compactness, we rewrite the Hamiltonian in terms 
of the Pauli matrices while rotating basis 7r/2 from the 
basis of localized states in the left and right dots to 
bonding-antibonding basis: 



p(t)e t/T0 =J2— dt k dtk-i ... dh 
k=o T ° ^° ^° Jo 

x f dW(X k ) f dW(X k ^)... 
Jx Jx 

dW{X )p{t 1 t k ), (8) 



TiA he HJt ,w > 
n = —o z + -a x + -^-tr z X(t) 



a x X(t). (4) 



x 



where dW{X) is the distribution of X with the probabil- 
ity of p u for X = 1 and that of pd for X = — 1, with the 
constraint p u + Pd = 1 • The density matrix before the 
ensemble average, p(t,tk), is given by 

p(t,t k ) = S{X k] t,t k )---S(X 1 ;t 2 ,t 1 )S(X 0] t 1 ,0) 

xp^S-^Xo; t 1 ,0)S- 1 (X 1 ;t 2 , h) ■ ■ ■ S^(X k ;t, t k ), (9) 



In the following, we consider only the case of a sym- 
metrical coupled-dot system, e = 0, in which the effect 
of the bias fluctuation due to the dephasing starts only 
from the second order, Jg, as shown in the following, 
and is less effective in the perturbation regime. As re- 
duced Hamiltonian Eq. (4) suggests, the present results 
can also be applied to other quantum two-level systems 
in which telegraph-type fluctuation exists (see Sec. VI). 



III. METHOD 

We are interested in the time-evolution of the qubit 's 
two- by-two density matrix, p(t), with an arbitrary ini- 
tial condition at t = 0, p(0). If BCF is absent, starting 
with the initial condition that the left dot is occupied, 
for example, the density matrix at t — is given by 



p(t = 0) 



(5) 



where W^(f ) is the matrix of rotation from the localized 
basis to the bonding-antibonding basis. At time t, 



p(t) = e 



1 




l+cos At 
2 

—i sin At 



i sin At 

l-cos At 



(6) 
(7) 



where the bases of the inner matrix are the left and right 
occupancy states. Therefore, the density matrix shows 
Rabi oscillation with frequency A. 

In the following, we keep the matrix indices in the 
bonding-antibonding basis. To examine the instanta- 
neous potential change, we use the method of stochastic 
differential equations [21]. The density matrix averaged 
over all possible sequences of telegraph noise can be rep- 
resented as a series: 



where S{X; t, t 1 ) — S{X 1 t — t') is the unitary time evolu- 
tion operator, which is determined by 



fo tS&'t) =H (X)S(X,t). 



The explicit form of S(X, t) is given by 



(10) 



S(X,t) 

_ ( cos \at - i{ ^ +J J X ) sin \at -i^r sin \at 

-i^sm\at " cos \at + { A +-^x sin i fflf J ■ 

(11) 



where a = \/(A + JtX) 2 + Jg. Equation (8) can be 
rewritten in terms of the integral equation 

p(t)e T l TQ = [ S(X;T,0)p(0)S- 1 (X;T,0)dW(X) 
Jx 

+ T j e ' /T ° / S(X;T 7 t)p(t)S- 1 (X;r,t)dW(X)dt. 

(12) 

Using S(X;t, t), we define matrix R lm (r,t) as follows: 

Rfk(r,t)= f S ik (X;T,t)S^(X;r,t)dW(X), (13) 
Jx 

B i ^{r,t) = {R k fr{r,t). (14) 



We can then reduce Eq. (12) to the following compact 
form: 

Pim(r) = e- T ^TT[R lm (T,0)p(0)} 

+ £ /; dt exp[-^]Tr[i?"»(r, t)p(t)\. (15) 
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IV. FLUCTUATION IN TUNNEL COUPLING 

First, we consider the case of fluctuation in tunnel cou- 
pling (Jt ^ 0, Jb = 0). Since the interaction Hamilto- 
nian commutes with Hgb, the environment leads to pure 
dephasing without energy dissipation. We derive S(X,t) 
and from Eqs. (11) and (13) as follows: 



S kl (X, T -t) = exp[-(A(r - t) - J T X(r - t))(-l) k ]6 kl , 



R?k m (r -t)= [ cxp[»A(r - t) + iJ t X{t - t)] 
x 5 nk S hn dW(X). 



(16) 



(17) 



As a result, the off-diagonal element of the density matrix 
obeys the following integral equation: 

<t(t) = e- r / To i?(T) + — ( R{t- t)exp[-?^](T(t)dt, 
to Jo r o 

(18) 

where cr(r) = e~ lAr /0i 2 (t) /I/O12 (0) is a normalized off- 
diagonal element of the qubit density matrix measured 
in a Rabi oscillation frame, and 

R(t - t) = p u exp[z J T (r - t)] + p d exp[-i J t (t - t)]. 

(19) 

Equation (18) can be rewritten as a differential equation: 
o"{t) + -<t'{t) + (J| + i{^2±)J T ) a {T) = 0. (20) 



TO 



The initial conditions are ct(0) = = pi 2 (0)/|pi 2 (0)| 
and cr'(0) =0. As a result of coupling with the en- 
vironment, the off-diagonal element of the density ma- 
trix decays as a function of time. When the real part 
of the two roots of the characteristic equation of Eq. 
(20) almost completely degenerates, the short-time be- 
havior for t < min(-j^, 3r ) (initial regime) is not a sim- 
ple exponential decay. In this initial regime, the off- 
diagonal element of the density matrix shows Gaussian 
decay, a(t) ~ <r(0)(l - J T t 2 /2 + ■ ■ •) ~ a(0)e- J ^ 2 ^, ir- 
respective of the asymmetric probabilities, p u and pd- 
The decay of the off-diagonal element of the density 
matrix becomes exponential for the asymptotic regime, 
t > (l/(^/|^ -4J| |), when p u = p d . For J t t < 1, 

this criterion is obtained when one of the two exponen- 
tial decay terms becomes negligibly small. For Jtt 3> 1, 
the time constants of the envelope of the two dumped 
oscillating terms are the same, l/(2ro). Exponential de- 
cay appears after the inverse of the oscillation frequency: 



l/y/jljf - ^|. When p u = p d and J t t = \, one 



ob- 



tains, a(t) = e~~(l + 2^r), where the dephasing can 
never be a simple exponential decay. 

The time constant of this exponential decay corre- 
sponds to the dephasing time, T 2 . For p u ^ p d [18], 



I to V T o 



4J| 



u{ Pu_jn, )Jt] (21) 

TO 




IJtI^o 

FIG. 2. Dependence of dephasing rate T 2 1 / 1 JV | on 
tunnel coupling constant |Jt|to for various values of p u 
and p d with Jt / and Jb = 0. 

Figure 2 shows the | Jt|to dependence of dephasing rate 
T^ 1 for p u = 0.5, 0.6, 0.7, 0.8, and 0.9. In the two limits 
of weak and strong coupling, we have 



(1 - (p u -pd) 2 )J T T , 
(1 - \p u -p d \)/2T , 



l/r > \Jt\- 
1/t < \Jt\- 



(22) 



Namely, for a fixed \Jt\, when ro increases from 0, T 2 _1 
first increases and then decreases. It has a single maxi- 
mum of (1 - y/\p u -pd\)\Jr \ when r = 1/(2| J T \); there- 
fore, for any parameters, T 2 > 2r . Changing weight p u 
to make it more asymmetric reduces T 2 ~ 1 . In the limit 
of p u — ► or 1, the environment is nearly frozen, so the 
dephasing time becomes infinity. 



V. FLUCTUATION IN ASYMMETRIC BIAS 

Next we examine the effect of bias fluctuation ( Jt = 
0, Js 7^ 0). We consider only the case of symmetrically 
weighted telegraph noise (p u = p d = 1/2) for simplicity. 
In this model, H q t and the interaction Hamiltonian do 
not commute and the dephasing process is accompanied 
by a relaxation of the population. The unitary operator 
is thus 



4 



S(X,t) 

cos \Vtt — i — sin ^flt — i^f^siniftt > 

-i^fpsin|m cos \Slt + i% sin \Q.t ' 1 ^ ' 



where Q, = \J A 2 + J|, is the nutation frequency. Af- 
ter averaging, the few non-vanishing elements of matrix 
Kjg{t) have the form 

R{l(T-t) = B%(T-t) = l-P 12 (T-t), 
R\\{t -t) = R 22 (t -t)= R\ 2 (t -t) = RI\(t - t) 
= Pi 2 (r-t), 

R\1{t ~ t) = R\\*{t -t) = [P 12 {t -t) + cosO(r - t) 
+ ^sinn(T-t)], (24) 



where 



Pi 2 (T-*) = ^fsin ; 



J% . 2 fi(T - t) 



(25) 



A. Analytic solutions 

Differential equation (27) with initial condition Eqs. 
(28 - 30) has explicit solutions: 



3 i 



x (dc cos <p + id s sin 0) + ie "sin^, (31) 
where u — t/t . For i = 1 — 3, the coefficients are given 

by 



a c = ^-((i + /? 2 )A 2 

+ (l + /3 2 -2a 2 )A l + (l + /3 2 ) 7 2 
C ls = -J-a(2A 2 + 3A 4 + l + /3 2 ), 



3a 2 



A,; 



(A 2 + 2A, + 1 + 7 2 )(3A 2 + 2A, + 7 2 ), 



(32) 
(33) 
(34) 



and the A's arc the three solutions of 



The time evolution of the diagonal element of this ma- 
trix is determined by the following differential equation 
[21,22]: 



2 

TO 



1 



n'"(r) + f n"(r) + (4 + A 2 + J 2 B )n'(r) + ^n(r) = 0. 



(26) 



7 2 

TO 



Here, n(r) = pn — p 22 , and the initial conditions are 
n"(0) = -Jfn(0) and n'(0) = 0; and n(0) depends on 
the initial population having an arbitrary value between 
1 and -1. The time constant of the exponential decay of 
n(r) is the relaxation time of the population, T\. When 
\J B \ < A,1/t (weak coupling), Tf 1 - J|r . When 
\J B \ > A, 1/to (strong coupling), Tf 1 - l/2r . 

The obtained differential equation for the off-diagonal 
element is 

<{r) + -< 2 (t) + (i + A 2 + J 2 + i-K 2 (r) 



T() 



TO 



A A 2 7 2 7 2 
+ + — + ^-)a 12 (r) = ^-a 21 (r), (27) 



to 2r ' 



2t 



where ct 12 (t) = Pi 2 (t)/|pi2(0)| is the normalized off- 
diagonal element of the qubit density matrix. The initial 
conditions are 



<7'/ 2 (0) = (-A 2 
a[ 2 {0) = -iAa 12 (0), 
<7i 2 (0) = e**, 



4fw 2 (o) 



-<7 2 l(0), 



(28) 

(29) 
(30) 



where </> is the initial phase of the off-diagonal density 
matrix element. 



A 3 + A 2 + 7 2 A + a 2 = 0, (35) 

where a = t A, (i = t Jb, and 7 2 = a 2 + f3 2 . 

We will show one special case and three asymptotes. 

(3 = isolated system: Since Eq. (35) has solutions 
Ai = — 1 and A 2j 3 = Ha, the coefficients arc determined 
as C\ c = C\s = 0, C 2c = C?, c = 2 and C 2s = C^ s = — { . 
Therefore, we get (J\ 2 (u) = e - lau + l< t> ^ which is simply a 
natural rotation of the off-diagonal element of the density 
matrix. 

a, (3 <C 1 fast modulation: In this asymptotic case, 
the roots of Eq. (35) are Ai ~ — 1 + 1 and A 2 3 

-h 



— \P 2 ± ia. After determining coefficients Cj, we obtain 



a 12 {u) 



~ e 



-u+i4 



(36) 



Therefore, the off-diagonal element of the qubit density 
matrix decays exponentially at a rate of T 2 _1 ~ | J b tq. 

2 

a,(3*^>l slow modulation: The roots arc Ai ~ — ^ 



and A 2: 3 



2T 2 



±7*, so 3ta 12 (u) 



r(/?V 



a 2 e COS7U) cos0 + 



-e ^ 



sin 7U sine 



3(Ti 2 (u) <~ e ^i jU {— — sin7ucos0 + cos7usin0) 
particular, for 1 <C a <C (3 (strong coupling limit), 

a 2 

ci 2 (t) ~ e ^ u cos0 

+ ie~^ u cos/3usin</>. 



and 
In 



(37) 



These apparent different time dependences between the 
real and imaginary parts stem from the choice of coupling 
in the form Jb<J x X/2. If we choose the form Jb<J v X/2 
instead, the time dependences of the real and imaginary 
parts are interchanged. In this strong coupling limit, the 
time evolution of 0i 2 (t) explicitly depends on its initial 
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phase, <fi, not the simply like e 1 ^ . Therefore, as will be 
discussed later, if there are several such impurities, the 
total time evolution of o-i 2 (t) is not the simple product 
of each impurity's contribution. For 1 <C (3 <C a (weak 
coupling), we have 



0.8 



(38) 



In this case, we have exponential decay with T 2 ~ 
Jl J_ 

For a, 1 » (3 weak coupling or preservative 



S 0.6 



regime, have Ai ~ —1 — a t +1 i, and A 2 



r)»- 



0- 

.2 

— £ -U-\-l4 



± 



0.4 



(! + 2(J + 1) K S0 CT 12 

Therefore, we again have exponential decay with T 2 ~ 

J B r a 



T =10 1 °S 

A=10 1 V 

j R =ioV 



exp(-J B t t/2) 

exp(-J B 2 T T) 
n(x) 



_i 1 i_ 







20000 



,, . . 2 -5s . This coincides with the Redfield result, which 
was obtained by perturbation theory and is justified in 
the weak coupling case, | Jb\ <C 1/t [23,24]. Taking the 
limit a <C 1 further, we restore the result for a, (3 <C 1, 
fast modulation. 

To summarize, the dephasing rate is given by 



x/x, 







FIG. 3. The t/to dependency of density matrix of n(r) 
and |<7i2|(t) (solid curve) when Jt — with Jb = 10 s s _1 
and A = 10 10 s _1 . Dotted lines are analytically obtained 
asymptotic curves, which is almost identical to the solid 
curves. 



{ J 2 B r /2, 



T^ 1 = { l/2r , 



4/2 A 2 r , 
[ Jir /2(1 + A 2 r 2 ) 



for A,\J B \ « l/r 
for the real part and 
l/r < A < |J B | 
for the imaginary par 
and 1/to < A < |J B | 
for 1/t <C \J b \ < A 
For|J B |«i,A. 



(39) 



In all regimes, T 2 > 2r . 



B. Numerical results 

Here, we show the results of solving Eqs. (26) and (27) 
numerically. Figure 3 shows the t/tq dependence of n(r) 
and |o"i2 (t) | when \Jb\ <C A, 1/ro (weak coupling) along 
with the asymptotic curves obtained analytically. 

Figure 4 shows the t/to dependence of n(r) and 
|<Ti2(r)| in the case of strong coupling. It also shows 
the asymptotic envelope for n(r). 

Figure 5 shows the t/t dependence of 5Rcti2(t) in the 
case of strong coupling. It also shows the asymptotic 
curve. As shown in Figs. 3, 4, and 5, in the two con- 
trasting limits, the numerical and analytical results co- 
incide very well. It should be noted that we do not find 
Gaussian decay of the off-diagonal element of the density 
matrix for the initial regime, in contrast to the fluctua- 
tion in the tunneling coupling constant. 
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exp(-2A 2 t/J B 2 x ) 



exp(-t/2x ) 
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t/T 



FIG. 4. The r/ro dependency of density matrix of n(r) 
and |<Ti2|(t) (solid curve) when J T — with J B = 10 12 s _1 
and A = 10 10 s _1 . Dotted line is analytically obtained 
asymptotic envelope curve. 
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exp(-A 2 x/J B 2 x ) 
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A=10 1 V 

j R =io 1 V 1 
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x/x, 



1000 







FIG. 5. The t/to dependence of real part of density 
matrix cti2(t) (solid line) when Jt = with Jb = 10 12 s _1 
and A = 10 10 s _1 . Dotted line is analytically obtained 
asymptotic curve. 
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10 ,u 10" 

^o(s) 

FIG. 6. The r dependence of T' 1 and 2TT 1 when 
Jt = with J B = 10 8 s _1 and A = 10 10 s _1 . 

Next, we examine the To dependence of the relaxation 
rate, particularly for \Jb\ < A. Figure 6 shows the To 
dependences of T-f 1 and T 2 _1 when A = 10 10 s _1 and 
Jb = 10 8 s _1 . The limits of the long and short To fit well 
with the analytical asymptotic given by Eq. (39). The 
tq dependence of the relaxation time is fitted as 



Tr 1 = 27T 1 = - 
l 2 1 



+ A 2 7 



(40) 



When To < A -1 , the relaxation rates increase with To- 
When To > A -1 , the rates decrease with an increasing 
To. The shape of the dephasing rate as a function of tq 



is explained as follows. When To < A -1 , many dephas- 
ing events occur during one Rabi oscillation cycle, each 
event leads to collective disturbance. Because the long 
time constant of telegraph noise leads to large fluctua- 
tions in the variance of the rotating angle in the Bloch 
sphere during Rabi oscillation, the dephasing time de- 
creases with an increasing To. When To > A -1 , Rabi 
oscillation occurs over more than one cycle in time To 
and, in this regime, each dephasing event is independent. 
Hence, the dephasing time increases with To. The maxi- 
mum around tq ~ A -1 is a kind of resonance. 



VI. DISCUSSION 

We considered the effect of electrostatic disturbance 
due to background charge fluctuations (BCFs). To sum- 
marize our findings, for pure dephasing {Jt ^ 0,Jb =0), 
2to < T2 and T\ = 00. For dephasing with relaxation of 
the population (J T = 0, J B + 0), 2t < T 2 = TF X . 

Next, we discuss the relationship between the exper- 
imentally observed T2 and our results. In the present 
study, the dephasing time with a single background 
charge was found to be longer than the time constant 
of the telegraph noise for both tunneling and bias fluctu- 
ations. The observed time constant of a dominant ran- 
dom telegraph noise is about 30 /is or longer [10,12], so 
a rather long dephasing time is expected. However, in 
another experiment, the dephasing time was about 1 ns 
[25]. Therefore, a single telegraph noise source may not 
be enough to explain the experimental results; we should 
thus consider the effect of many impurity sites or other 
additional effects. 

We consider the effect of many impurities for the case 
in which there is no correlation between background 
charges. Flucuation in the tunnel coupling arises from 
the modulation of the wave function in the coupled dots. 
The gradient of the electrostatic potential around the 
tunneling barrier, which comes from an electron or hole 
located at an impurity site, leads to a change in the tun- 
nel coupling as well as a change in the asymmetric bias 
[26]. For a charge state that couples with a coupled- 
dot system symmetrically, the pure dephasing event is 
critical. For dephasing accompanied with relaxation of 
the population, the unitary operators of each impurity, 
which lead to dephasing of the qubit, are not commut- 
ing. However, by neglecting the higher order JgTo's in 
the dephasing rate, we can take the ensemble sum of the 
effect of each charge state. In such a weak coupling case, 
(JbTq -C 1), we can use Eq. (40), and the simple sum- 
mation of the dephasing rate is expressed by [20] 

T£ = J dJ B Pz(J B ) J dWP{W)T^{J B ,W) 
_ k B Tir < J% > s 



W 



4A 



(41) 
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where Wo is the distribution width of the thermal acti- 
vation energy of the charge states, Ps(Jb) is the dis- 
tribution function of J B , which depends on the rela- 
tive position between the qubit and impurity site, and 
< J B >s (= J cLJbP?,{ Jb) J%) i s the sum over the ran- 
dom impurities. In the second equation, we assume 
uniform distribution of the activation energies of the 
background charges, P(W) — 1/Wo; for typical cases, 
Wq/UbT is approximately 23 [20]. Use of the perturba- 
tion method showed that the dephasing rate of a Joseph- 
son charge qubit in terms of JbTq is proportional to the 
inverse of Ej in the limit of Eq = [14,15], where Ej 
and Ec are the Josephson coupling constant and charg- 
ing energy, respectively. This is similar to the estimate 
of T 2 ^ in Eq. (41), where Ej/Ti corresponds to A. With 
a larger A and lower temperature, the dephasing rate is 
lower. 

We next estimate the magnitude of the fluctuations, 
Jb- The asymmetric bias fluctuation comes from asym- 
metric coupling between the two dots and the back- 
ground charge, which is in the form of a dipole inter- 
action, Jb oc e 2 dcos8/r 2 for rf<r, where d is the dis- 
tance between the two dots, r is the distance between 
the coupled-dot system and the background charge, and 
9 is the angle between them [27]. Therefore, for a smaller 
qubit or a charge located far from the qubit, the effect of 
bias fluctuation should be less important. The dephasing 
rate is proportional to < J B >£, which is estimated as 



4ne r e hr 



cosO,) 2 ~ ( 



e 2 d 



4:ne r e h Jr m (d) 

. „cos 2 9 
x / sin 9 ^d92n 



d , 9 4tt Nj 



r dr 



= (t—t) 



4e r e h 3 r m (d) 



(42) 



for the impurity sites where the dipole approximation is 
appropriate, where r m (d) is the radius beyond which the 
dipole approximation is valid, which depends on d, Ni 
is the density of impurity sites, and e r is the relative di- 
electric constant. Therefore, the total dephasing time is 
well defined. The quality factor of a quantum logic gate 
is defined by the ratio of A to : 



Q 



w 



4A 2 



-2S 



k B Tir < J 2 > s ' 



(43) 



which represents how many gate operations can be done 
before the quantum coherence vanishes. From Eq. (43), 

A 2 

we conclude that a large , 2 is needed for quantum 
computation. We estimated Q using d — 0.3 fim, r, m = 1 
/im, A = 200 fieV (characteristic parameter values for an 
experiment in which the quantum mechanical coupling of 
the dots was observed in the frequency domain [5] ) , and 
e r = 12.5 (for GaAs). To enable quantum error correc- 
tion, the lower bound of the necessary gate quality factor 



was roughly estimated as Q > 10 6 [28]. Thus, density of 
charge states should be less than 5 x 10 6 cm~ 3 for fab- 
rication. If there is a correlation between impurities (a 
screening effect), dephasing will be suppressed in general 
[29]. It should be noted that there might be strong de- 
phasing from the nearby impurities for which dipole ap- 
proximation is not adequate, even if these arc only a few 
impurities, (a few in this case). The non-commutativity 
between the qubit Hamiltonian and environment Hamil- 
tonian and the qubit backaction make it difficult to eval- 
uate the dephasing rate for strongly coupled background 
charge fluctuations in the asymmetric bias case [14]. 

Finally, we discuss the Josephson charge qubit sys- 
tem [1]. Under an appropriate condition (single-electron 
charging energy Ec much larger than Josephson coupling 
energy Ej and temperature fcgT < Ej ) only two charge 
states in the Cooper pair box (CPB) are important, and 
the Hamiltonian is given by 

rr Ej SE C h^JaXj 

H = —<r x + —<?z + 2 CTz ' ( U > 

where SEc = 4Ec{Qt/c- — 1) is the energy difference 
between the two charge states, and Qt is the total gate- 
induced charge in the box. The two-charge-state basis is 
expressed using Pauli matrices, and HJc is the coupling 
strength between the qubit and the background charge, 
which induces fluctuation in the charging energy. The 
Ej/h corresponds to the asymmetric bias e, and Ec/h 
corresponds to A. Here, Ec — 122 /ieV, and Ej ~ 34 
/j,eV [16]; if we can neglect Ej, the pure dephasing event 
is critical. In pure dephasing, the effect of a large num- 
ber of impurities is obtained by simply summing the de- 
phasing rates, because H q b and interaction Hamiltonian 
commute. When the background charge and CPB in- 
teract, the charging energy in the CPB fluctuates. The 
spectrum of the fluctuation is given by 

Sae(co) = J dte wt Y,K 2 < JaXiWJcjXjiO) > 



= J dJcP( Jc) J droP(r )^^ 



k B T nh 2 <J C > 



Wo 



2u) 



(45) 



where < J c >= f dJcP(Jc)Jci an d we take an en- 
semble average over the activation energy, as was done 
in the coupled-dot system. The spectrum of the charg- 
ing energy fluctuation was experimentally found to be 
S A e{uj) = (^f 2 ) 2 ^ whereas (1.3xl0- 3 e) 2 [16]. From 
this estimation, < J c >~ 4.6 x 10 23 s~ 2 for 20 mK. For 
an initial regime, the envelope of Rabi oscillation shows 
Gaussian decay, namely, the off-diagonal element of the 
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density matrix is given by pi2{t) — Pi2(0)exp(— < > 
t 2 /2)e~ lEjt / h . The rate of Gaussian decay is given by 
\f< Jq > /2. For the above value of < Jq >, the time 
constant of the Gaussian decay is given by 3.6 x 10 2 ps, 
which is consistent with the experimental finding of 150 
ps 7 [16]. Note that in the Gaussian regime, the time con- 
stant does not depend on the temperature. Numerical 
calculation [14] also suggests this type of Gaussian decay 
for the pure dephasing case. 

At the charge degeneracy point, namely 5Eq = 0, de- 
phasing with relaxation of the population occurs. In this 
case, the dephasing rate is estimated using the same value 
of < j£ >: 



! = 7T k B T <Jji> 
2 4 Wo Ej • 



(46) 



so T2 ~ 0.28 ps for T = 20 mK. In a recent experiment, 
a longer coherence time of 0.50 /is was found when the 
saddle point of the ground state energy was used as a 
function of Q t and the flux [30]. 



VII. SUMMARY 



We examined the effect of the fluctuation of a single 
charge in an impurity site on a qubit. Using the method 
of stochastic differential equations, we calculated the 
time evolution of the ensemble averaged density matrix of 
the qubit and obtained analytical results for various con- 
ditions. The dephasing time, T 2 , was always longer than 
the time constant of the random telegraph noise for both 
tunneling and bias fluctuations. For bias fluctuation, T2 
was twice the relaxation time of the population in the 
weak coupling case. To suppress the bias fluctuation, 
the coupled dots should be positioned closer together or 
the tunnel coupling should be made stronger. We also 
investigated the case in which many impurity sites are 
distributed and examined the gate quality factor. For 
pure dephasing, which corresponds to a Josephson charge 
qubit experiment, the Gaussian decay of the off-diagonal 
element of the density matrix dominated. The present re- 
sults can be applied to other quantum two-level systems 
in which there is telegraph-type fluctuations. 
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